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Abstract 

The Inozemtsev Hamiltonian is an elliptic generalization of the differential oper- 
ator defining the BCjv trigonometric quantum Calogero-Sutherland model, and its 
eigenvalue equation is a natural many- variable generalization of the Heun differential 
equation. We present kernel functions for Inozemtsev Hamiltonians and Chalykh- 
Feigin-Veselov-Sergeev-type deformations thereof. Our main result is a solution of a 
heat-type equation for a generalized Inozemtsev Hamiltonian which is the source for 
all these kernel functions. Applications are given, including a derivation of simple 
exact eigenfunctions and eigenvalues for the Inozemtsev Hamiltonian. 
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1 Introduction 



Integrable models in quantum mechanics are closely related to the mathematical theory of 
special functions. A famous example are Calogero-Moser-Sutherland models which describe 
an arbitrary number of identical particles moving in one dimension under the influence of 
particular one- and two-body potentials [H [2j [3j 0]. The Hamiltonians of these models 
are differential operators that have eigenfunctions providing natural many-body generaliza- 
tions of the classical orthogonal polynomials. For example, the so-called BCat trigonometric 
Calogero-Sutherland model has energy eigenfunctions given by many-variable Jacobi poly- 
nomials; see e.g. [HE]. 

* langmann@kth.se 
ttakemura@math.chuo-u.ac.jp 
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The BCjv trigonometric Calogero-Sutherland model has an elliptic generalization defined 
by the Hamiltonian 



j=l 3 u=0 (!) 

+ ^ 2 A(A ~ l ) {P( X J ~ X k) + p( x j + x k)} 
±<3<k<N 

with p(x) the usual Weierstrass elliptic function with periods 2ui and 2co> 3 and where we 
use the notation 

u = 0, u 2 = -wi - uj 3 , 

here and in the following^] This Hamiltonian depends on the particle number N and 
five coupling parameters g , gx, g 2 , #3, A. It defines the natural quantum-analogue of a 
classical Liouville integrable system first presented by Inozemtsev in [7] , and we refer to this 
quantum-many body system as the Inozemtsev model. The integrability of the Inozemtsev 
model was partially established by van Diejen [8], and Oshima j9] described its commuting 
operators (higher-order Hamiltonians) completely. In the following we sometimes write 
Hn (x; {g u }l =0 , X) for the Inozemtsev Hamiltonian in (jTJ), to indicate the argument x = 
(xi, . . . ,2at) and the coupling parameters. Note that this differential operator allows for 
a quantum mechanical interpretation only if one assumes that uj% > 0, — 1U3 > 0, and 
that all coupling parameters are real. However, many of our results hold true with lesser 
restrictions. 

The Inozemtsev model is interesting already in the one-variable case: the eigenvalue 
equation Hiip(x) = Eip(x) of the differential operator 

Q2 3 

Hi = -— + ^2g v (g v -l)p(x + u v ) (2) 

is equivalent to the Heun differential equation, which is a standard form of a second-order 
Fuchsian differential equation with four singularities and a topic of current research in 
special function theory; see [101 EJ EG] . The Heun differential equation, and the differen- 
tial equations of its confluent type, appear in several physics contexts, including quantum 
mechanics, general relativity, models of crystal imperfections [TT], and the AdS/CFT corre- 
spondence p3] . In this paper we present generalizations of two known functional identities 
involving the Heun differential operator in (j2j) and functions that are products of powers 
of Jacobi theta functions 9 u+ i(x), v — 0, 1,2,3 (see Appendix IA.1I for precise definitions). 
The first known identity is as follows: The function 

3 

* 1 (x) = He u+1 (x) 9 » (3) 

obeys the equation 

2(00 +9l + ga + ga )- + H 1 -C 1 V ^i(x) = (4) 



1 Our notation for elliptic functions is as by Whittaker and Watson [14], except that 0J2 is denoted by us 
as W3; for the convenience of the reader we collect the definitions of the functions we use in Appendix lA.il 
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with a known constant C\ (see (126|) ). where /3 = 2U1U3/ (iri) {uj\ is fixed). This non- 
stationary Heun equation appears in several physics contexts, including the Wess-Zumino- 
Witten model [15], the eight- vertex model [16] and Liouville field theory [17] . The second 
known identity provides a kernel function for a pair of Heun differential operators (recall 
that a function F(x, y) of two variables x and y is called a kernel function of two differential 
operators D(x) and D(y) if [-D(x) — D(y) — c]F(x, y) = for some constant c): It is known 
that the function 

is a kernel function of the Heun differential operators H(x; {g u }l=o) an d ^(y> {9v}l=o) 
provided that g v = A — g u and A = (g + g\ + g 2 + g?)l2 [HI 120]. In fact, we find that, for 
arbitrary A, the function in ([3]) is a generalized kernel function of these differential operators 
in the following sense, 

d 1 

2(0o + 9i + 92 + 93 ~ 2A)— + H^x; - H^y; {g„}l =0 ) + C M j * 1A (x, y) = 

(6) 

with a known constant C^i (see (I30j) ). Moreover, our results include another generalized 
kernel function for two Heun-type differential operators: the function 

3 

^(x, y) = 6 1 (x- y)B x {x + y)\{ v+1 {x)*'0 v + 1 {vfi> (7) 

obeys 

2(g + gi + g 2 + g3 + 2)— + H x {x- {^}Lo) + AfliG/; {gX=o) + Ci.i | *i,i(*> 2/) = 

(8) 

provided that ^ = (2^ + 1 — A)/ (2A), for some known constant Ci,i (see ( 134"1) ) and arbitrary 
A / (I. 

In this paper we present and prove many-variable generalizations of the identities in the 
previous paragraph: we obtain a generalizations of the non-stationary Heun equation in 
©-(BJ, an d of the two kinds of generalized kernel function identities in (J5])- (J6j) and ©-([8]), 
to Inozemtsev Hamiltonians, for arbitrary particle numbers; see Corollaries 13 . 1\ 13721 and [3731 
respectively. Note that, in the latter two cases, the two Inozemtsev Hamiltonians can have 
different particle numbers N and M. The most general kernel function identity we obtain 
is for a pair of differential operators 

H N)R {x, x- {g„}l =0 , A) =H N (x; {g v }l=^ A) - \H^(x; {(A + 1 - 2g v ) / {2X)} 3 V=0 , 1/A) 

N N ^ 

+ ^2^2(1 - A) {p(xj - x k ) + p(xj + x k )} 

j=l k=l 

and H M ^(y,y; {A — g v }l =1 , A), with Hn(x; {g u }l =0 , A) in p]), for arbitrary particle numbers 

N, N, M, M; see Corollary 13.41 All our results are special cases of a many-variable 
generalization of the non-stationary Heun equation in (E]) (J4]) to a generalized Inozemtsev- 
type differential operator where all particle mass parameters can be different and where the 
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interaction strengths depend on these "masses" in a particular way; see Theorem 12.11 We 
refer to this as source identity since it is the source of all our other results: the latter are 
obtained in a simple way as special cases. Remarkably, direct proofs of these special cases 
are often more complicated than the proof of the source identity. 

It is known that integrable quantum mechanical models of Calogero-Moser-Sutherland 
type allow for deformations that share many of their beautiful mathematical properties 
[221 1231 124"] . and for the BC^ trigonometric Calogero-Sutherland system this deformation 
corresponds to the trigonometric limit of the differential operator in ([9]). It thus is natural 
to conjecture that the generalization of this deformation to the Inozemtsev model is given 
by ©. 

The differential operator defining the Sutherland model (= A^_i trigonometric Calogero- 
Sutherland model) has a well-known kernel functions which can be used to construct the 
eigenf unctions and eigenvalues of the Sutherland model [25], and this approach can be 
generalized to the elliptic case [26] and to all quantum Calogero-Moser-Sutherland models 
associated with classical orthogonal polynomials [BJ. Moreover, the latter results allow for a 
natural generalization to the deformed models [BJ. Our results in the present paper provide 
the starting point to generalize these results to the Inozemtsev mode and its deformation 
in ©. 

The first example of a source identity was found for the Sutherland model by Sen [27J , 
and that this identities can be used to obtain kernel functions for the Sutherland model and 
its deformations was pointed out in [28] . Source identities for all quantum Calogero-Moser- 
Sutherland models where obtained and used to derive kernel functions in |6J. A source 
identity allowing to derive kernel functions for the elliptic generalizations of the Sutherland 
model and their deformations was presented in [25] • The present paper generalizes, to the 
elliptic case, results previously obtained in [BJ. 

We mention four further topics for future research suggested by the results in the present 
paper. First, kernel functions of Calogero-Sutherland models can be regarded as a natural 
quantum analogue of Backlund transformations found by Wojciechowski [30]; see also [3Tj . 
This suggests that our results can provide Backlund transformations for the classical version 
of the Inozemtsev model. Second, kernel functions for the Sutherland model have been used 
to construct Q-operators that allow to derive integral representations for the eigenf unctions 
of this model [52]. This suggests that it is possible to extend Sklyanin's separation-of- 
variable approach [33J to the Inozemtsev model using our results. Third, as pointed in [34J 
(see also [35]), some special cases of the Inozemtsev model are quasi-exactly solvable in 
the sense that the computation of a finite number, m, of eigenfunctions can be reduced to 
diagonalizing a m x m matrix (note that the Hamiltonian in Eqs. (2)-(4) in |34j is identical 
with the Inozemtsev Hamiltonian as in ([T]) for A = a, g = gi = g 3 = b, g 2 = —3b — 2[m + 
X(N — 1)]). We find that these cases are special with regard to our kernel functions (since, 
by applying the result in Corollary 13.31 to the parameters above and (N,M) = (N,m), we 
obtain A^,m = 0). This suggests that our results might shed new light, and possibly allow to 
extend, the results in [31]. Fourth, while many kernel functions for Ruijsenaars' relativistic 
generalizations of Calogero-Sutherland-type systems [3BJ are known [37], the results in [BJ 
and the present paper suggest that many more such relativistic kernel functions should 
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exist. We believe that a good strategy to find all such kernel identities would be to find 
relativistic generalizations of the source identities obtained in [6j [29] and the present paper. 



2 Main result 

As discussed, our main result is a heat-type equation for a Schrodinger-type differential 
operator. Our notations for elliptic functions is defined in Appendix lAl 

Theorem 2.1. (Source Identity): For M a positive integer, X, d u (v = 0,1,2,3,) and 
mj 7^ (J = 1, 2, . . . ,Af) complex constants, and Xj (J = 1, 2, . . . ,Af) complex variables, 
let 



= ( n n ^+i{xj)^A n e^xj - x K r^%{xj + x Ki 

\j=lv=0 J 1<J<K<M 



(10) 



and 



N 



n = J2— "ax 71 + J29-A9u,j - i)p(Xj + u v ) 

J=l J \ J u=0 / ( | | ; 

+ U,k{p(Xj-X k ) + P (Xj + X k )} 



kj<k<M 



with 



Then 



with 



1j,k = A(mj + m K )(\m,jm K - 1) (12) 
9u,j = mjd u + ^m 2 j. (13) 

{(4A|m| + 2\d\)— +H- £ )%(X) = (14) 

S =(2X\m\ + \d\) {Af - X(\m\ 2 + \m 2 \) - \m\\d\} ^ 

^i (15) 

+ \m\{(d d 1 + d 2 d 3 )e± + (d d 2 + d 1 d 3 )e 2 + (d d 3 + did 2 )e 3 } 

3 M M 

\d\ = d u , \m\ = ^irij, \m 2 \ = s ^^m 2 J . (16) 

u=0 J=l J=l 

Proof. Consider the differential operator 

k = Y,^tQjQj ( 17 ) 



with 



Qj - ± dX~ J +Vj ' Vj ~MX) dXj • 
5 



Using identities of elliptic functions collected and proved in Appendix IA.2I we find, by 
straightforward computations (details are given in Appendix IB]) . 

U = (4|m|A + 2M)-J- J-$ + n - So- (19) 

By definition, Qj§o = for all J. Thus "H$o — 0, and (|T9l) implies our result in (fT4")) . □ 

Under suitable restrictions on parameters (see below), the differential operator ~H in 
( fill has a natural physical interpretation as Hamiltonian describing Af distinguishable 
quantum particles with interactions, and the results above provide the exact groundstate 
and groundstate energy of this Hamiltonian. Namely, if cu\ > 0, — iw 3 > 0, A > 0, mj > 0, 
d Q > —\mj/2 and di > —\mj/2 for J = 1, 2, . . . , Af, d u real for v — 2, 3, and 4|m|A + 2|d| = 
0, then the Hamiltonian in (fTT|) defines a unique self-adjoint operator on the Hilbert space 
L 2 ([0, wi]-^) which has &o(X) as groundstate and So as groundstate energy. (This is true 
because, under these conditions, So and the potential terms in the Hamiltonian 1-i are real, 
the function $o(AT) is square- integrable, and the Hilbert space adjoint of Qj is equal to 
the closure of Qj. This, f fT7|) . f lT9|) . and the vanishing of <9/<9/3-term in f jT9|) imply that 
Yl j(QjV Qj I m J + So defines such a selfadjoint extension 7-L; see e.g. [38])- In the rest of 
this paper the self-adjointness of Inozemtsev-type differential operators will play no role. 



3 Special cases 

To state important special cases of our main result we use the following notation 

N ( d 2 3 \ 
j=i x i u=o (20) 
+ ^2 - !) {p( x j - x k) + p{xj + Xk)} , 

l<j<k<N 

* N (x;{g v }l=oA)= (nn^+ 1 ^')"j II O^-x^e^+x^ (21) 

\j=l u =0 ) l<j<k<N 

for x — (xx, . . . , Xjst) and complex variables Xj. We also use the abbreviations 

CO = {(0001 + 0203)61 + (0002 + 0103)e 2 + (0003 + 0i02)e 3 }, (22) 
101= 00 + 01 +02 +03- (23) 

We first state the many-variable generalization of the non-stationary Heun equation in 

Corollary 3.1. For N a positive integer, g v (v — 0, 1, 2, 3) and A complex parameters, the 
following holds true 

A N ^ + H N (x; {g u } 3 u=0 , A) - Cjvj ^n(x) = (24) 



with 

A N = AX(N - 1) + 2\g\ (25) 
Cn = 4^[1 - X(N - 1) - \g\]^ + Nc . (26) 

Proof. Set M = N, d v = g v - A/2 {y = 0, 1, 2, 3), (mj, Xj) = (1, Xj) for J = 1, 2, . . . , N 
in Theorem 12.11 and rename "H, $oPO, £o to Hn(x), ^^(x), Cn, respectively. Recall 
ei + e 2 + e 3 = 0, which implies {(d d 1 + d 2 d 3 )e 1 + (rf ^2 + e£id 3 )e 2 + (d d 3 + did 2 )e 3 } = c . □ 

The many-variable generalization of the generalized kernel function identity in (J5j) — (J6j) 
is as follows. 

Corollary 3.2. For N, M non-negative integers such that N + M > 0, g v (v = 0, 1, 2, 3) 

and A complex parameters, let g v = A — g v and 

9 NM (x,y) = ^y^Aj^^a^o.^ ^ (27) 

Then 

|^tv,a/^ + H N (x; {g v }l =0 , A) - H M (y; {g u } 3 u=0 , A) - C N>M \ ^ N<M (x,y) = (28) 

with 

A N>M = 4\{N-M-l) + 2\g\ (29) 

C NM = ^{(iV + M)(l - A) - (N - M)[(N -M- 2)A + \g\}}^- + (N - M)c . (30) 



Proof. Similarly as above, but now set TV = N + M, d v = g v — A/2 [v = 0, 1, 2, 3), and 
(mj,Xj) = 

in Theorem 12. 1L □ 



j = i,...,iv 

;-1,2/j-at), J = N+1,...,N + M 



The many-variable generalization of the generalized kernel function identity in (JTJ) — (JH1) 
is as follows. 

Corollary 3.3. For N, M non-negative integers such that N + M > 0, g v (v = 0, 1, 2, 3) 

and A 7^ complex parameters, let g' u = (2g u + 1 — A)/(2A) and 

(N M \ 
II II 01 ( x i - Vk)H*j + S/*) {9u}l=o, A)* M (y; {^}"=o> 1/A)- 

3=1 fc=l / 

(31) 

T/ien 

^v,m^ + H N (x; {g„}l =0 , A) + XH M (y; {g'„}l = o, 1/A) - CW.m }Vn,m( x >v) = ( 32 ) 



with 

A JV)M = 4A(iV-l)+4M + 2|£| (33) 

C NM = ^^{N+M-(N+M/X)[(N-2)X+M+\g\]-NX-M/X}—+(N+M/X)c G . (34) 
2 oj\ 



Proof. Similarly as above, but now set Af = N + M, d v — g v — A/2 {y — 0, 1, 2, 3), and 
(mj,Xj) = 

in Theorem 12.11 □ 



(l,xj), J=l,...,iV 
(l/\,yj- N ), J = N+1,...,N + M 



We finally state the generalized kernel function identity for deformed Inozemtsev Hamil- 
tonians. Note that all previous results stated in this section are special cases of this. 

Corollary 3.4. For N, N, M, M non-negative integers such that N + N + M + M>0, 
d v (v — 0, 1, 2, 3) and X ^ complex parameters, let g v — d v + A/2 (v — 0, 1, 2, 3) and 



H ( ^l(x, x) =H N (x; {A/2 ± d u } 3 u=0 , A) - XH^(x; {(1/2 =p d„)/A}jUo, 1/A) 

Af N 

+ 5^ 2(1 - A) {p(x 3 - - x k ) + p{xj + Xk)}, 

j=l fc=l 



(35) 



^(±)_ , ~s = *n(x; {A/2 ± d„}l^ A)^(x; {(1/2 =f 4)/A};Lo, 1/A) ^ 

x TT (n gl ^ 3 ' ~ r ^ m TT ~ r ^ ' ^ 

H vM nf=i ^ - 1 M nf =1 - r^v* 



Then 



A - 

^N,N,M,M go 

with 



,mqq + h n)<( x > ^ - H Mji(y> y) - c n,n,mm) ^n,n,mM x > 5 > y>y) = ( 38 ) 

^Vjv,a/,m = 4A(iV — M — 1) — A(N — M) + 2|p| (39) 
C n ,m M m = An >* 2 MM {N + N + M + M - \m\ [(\m\ - 2) A + \g\] - \m 2 \x] * + | m |c , 



2 I y '" u 1 J Wl 

|m| = N — M — {N — M)/X, |m 2 | = iV + M+ (iV + M)/A 2 . 



(40) 
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Proof. Similarly as above, but now set Af = N + N + M + M and 



{mj,Xj 



(l,xj), J=1,...,N 

(-1/A,£j_jv), J = N + 1,...N + N 

{-l,yj_ N _„), J = N + N + 1,...,N + N + M 

{V\Vj-n-n-m)> J = N + N + M + 1 j ...N + N + M + M 



in Theorem 12. 11 □ 

Note that H^fi(x,x) in ( 135]) is equal to H N ^(x,x; {g u }l =0 , A) in (jUJ), and H^^x, x) in 

(j3"5"|) is equal to H N j^(x, x\ {X — g u }l=o, A). Thus the general kernel identity in Corollary 13.41 
is a natural generalization of the one in Corollary 13.21 

We also note that the generalized kernel function identity on Corollary 13.41 is invariant 
under the following transformations, 

(N, N, M, M, {g u } 3 u=0 , A) (M, M, N, N, {A - g v } 3 u=0 , A) (41) 

(the variable names should also be changed correspondingly, of course) and 

(N, N, M, M, {g u }l =0 , A) (N, N,M,M,{(\+ 1 — 2 9v )/(2X)}l =0 , 1/A), (42) 

and these symmetries provide non-trivial checks of our computations. To be more specific: 
under the transformation in (j4"Tl) . the constants in (139]) and (j4"0]) change as ^riVMAf ~~ 
— -Ajv jv" m m> n m m ~* n m Mi an d under the transformations in (14"2"]1 they change as 
^n,n,m,m ~ A n,n,m,m/\ c n,n,m,m ~> ~ c n,n,m,m/^ consistent with the transformation 
properties of the r.h.s. of the generalized kernel function identity in (1381) . Note that not only 
(jHJ) but also (HSl) is a duality transformation (i.e. applying each of these transformations 
twice gives the identity). 



4 Applications 

Using the kernel functions obtained in the previous section it is possible to extend methods 
developed in 0,[26l|19] (e.g.) to construct eigenfunctions and eigenvalues of Inozemtsev-type 
differential operators. This section describes a general strategy and two simple examples. 
More systematic studies are left to future work. 



4.1 Integral transformations 

We explain how the kernel functions obtained in the previous section can be used to con- 
struct integral transformations that map a known generalized eigenfunction of a Inozemtsev- 
type differential operator to a generalized eigenfunction of another such operator. 

Let V Nt N,M,M = ^n,n,m,m( x ^iViV) be a generalized kernel function, H^l(x,x) and 
H^~'^(y,y) Inozemtsev-type differential operators, and A N ^ M ^ and C N fj MM constants 
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as in Corollary 13.41 If f(y,y) is a generalized eigenfunctions of the differential operator 

M,M V 



H[ \- (y,y) in the following sense, 



dp ' " !/ - w 
for some constant E, then ( 1381) implies 



a n,nmmjL + H lUy, v)~ E ) f(y, v) = ( 43 ) 



f{y,y) 



= f{y,y){H { M ^(y,y)^ Ni N tMt M} - ^^m.mI^m^^^)/^^)} 

^ d f d \ 

= ~J2dy- (f(y> y)-Q^N,N,M,M - ^N,N,M,Mg^f(y^ V)j 



(44) 



M 

dm,- 

k=l yK 

Integrating this with respect to the variables (y, y) over a suitable region C, one finds that 



f(x, x) = j V Nj n iM: m{x, x, y, y)f{y, y)d M yd M y (45) 
is a generalized eigenfunctions of the differential operator H^^(x,x), i.e., 

(^ A n,n,m,m-^ + H nJv( x > x)-E- C NflMj ^j f(x, x) = 0. (46) 

Note that a key point in the derivation of this result is that the region C is suitable in the 
following sense: First, the integral in (145]) has to be well-defined, and second, the integral 
over the total derivative terms in the last two lines of ( I44p must vanish (in general, Stokes' 
theorem implies that the latter is equal to an integral over the boundary of C). 



4.2 Example 1 

To be specific we assume throughout this section that u\ > 0, — ico> 3 > 0, and that Xj and Xj 
are real variables. As discussed, in this case the Inozemtsev Hamiltonian can be interpreted 
as a quantum mechanical model of a many-particle system. 



We consider the result in the previous section in the case M = 1, M = 0, and g u e {0, 1} 

M,M V 



(v = 0,1,2,3). Then H\ ' \~ (y,y) = -d 2 /dy 2 , and it is trivial to find solutions of (g3 



f(y) = exp(-ipy) and E = p 2 , with p an arbitrary constant. Moreover, 



$ Nf j,Jx,x,y) = ¥+Ux,x) I I I 9 1+ Jy) § " I J — J/ ^~ J ' Zl , (47) 
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A suitable integration region in this case is any path C in the complex y-plane such that 
n 1 o( x > ^> v)f(y) is analytic in some neighborhood C and such that 



f d ( d _ d \ 



(4* 



To find such a path we use that the Theta- functions 9i +u (x) can be expressed in terms of 
meromorphic functions 8i +u (z) of the variable z = exp(—2ix/R) as follows 



6 2 (x) = e T ^e^ R 9 2 (z), 9 2 {z) = ^° =o e™^n+i)( z -n + 

e 3 (x) = e 3 (z), e 3 (z) = i + EZie Tnin2 (z~ n + zn ) 



(49) 



(this is a simple consequence of the definition of the Theta-functions in (|66|) ). For p = 
(2n + go + gi)/R and n an integer, one thus finds that ^ N $ 1 (x, x, y) exp(— ipy) is a 
holomorphic function in the variable £ = exp(— 2iy/R) in the region 1 < |£| < q~ 2 , and 
that (14"5I) is satisfied for C the straight line from y = ie to 7r_R + ie, with e > such that 
exp(2e/i?) < g -2 (since this corresponds to a closed path in the complex £-plane where the 
integrand is holomorphic). Taking the limit e — > + we obtain the following. 

Proposition 4.1. For N, N non-negative integers such that N + N > 0, A a non-zero con- 
stant, g u G {0, 1} (v = 0, 1, 2, 3), and n an arbitrary integer, let d v = A/2 — g u , H^^(x, x) 

as in [35]) , and ^/^^(x,x) as in [36]) . Then the function 



f n {x, x) = *fy(x, x) lim / JJ 6i+ v (y) 

E_> J' le \u=0 



I&Mxj-yWifa + y)* 



(50) 



zs well-defined and obeys the equation 

A N,NX0§p + H n.U X > ^ ~ ( ^ + ~ 90 + ~ 9l)2/R2 ~ C W,w) f»( X > £ ) = ( 51 ) 

with 

A N ,N, lj0 = 4:\N-iN-2\g\, (52) 

A, 
"2 

Co = {(Mi + 92h)ei + (g h + 9ih)^2 + (Ms + 9xh)^}, 
\~9\ = ~go + ~9i + ~g2+h, \m\ = N-l-N/\, \m 2 \ = N + 1 + N/X 2 



C n ,n X0 = 1Z tt M {N + N+l- \m\[(\m\ + 2)A - \~g\] - \m 2 \X \ + \m\c , 



(53) 
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A noteworthy special case is N = 1, N = 0, and A = \g\/2. Then the function in (!50|) is 



with g v = \g\/2 — g u , and (l5"Tj) is the Heun differential equation 



(54) 



<9 2 3 

■q-2 + ^9u(g v - l)f>(x + u v ) - (2n + g + gi) 2 /R 2 ) f n (x) = 0. (55) 

It is worth noting that the integral factor on the r.h.s. in ( )50l) is proportional to the 
function f n (z,z), Zj = exp(—2ixj/R) and Zj = exp(—2ixj/R), defined by the following 
generating function, 

TTWO"" v = y fn{z,z)C n , (56) 



where the series on the r.h.s. is absolutely convergent in the region 1 < |£| < q 



-2 



4.3 Example 2 

We show how kernel functions can be used to transform Bethe ansatz solutions of the Heun 
equation obtained in [19] to eigenf unctions of Inozemtsev-type differential operators with 
arbitrary particle numbers N, N (this is a generalization of a result obtained in [12] for the 
case N = 1, N = 0). 

As shown in [19], for arbitrary non-negative integers n v [y = 0,1,2,3) and arbitrary 
constants E, the Heun differential equation 







- 3 _\ 
+ ^2n u (n u + l)p(y + u v ) — E 



dy 

has a non-zero solution which can be written as 



f(y) = (57) 



no+ni+n2+ri3 

exp(rey) | ( O^y + tj) 

j=X 

f{y) = 1 (y) n «9 2 (y) n '9 3 (y) n z6 4 (y) n 3 ' (58) 

for some constants tj (j = 1, . . . , n + ri\ + n 2 + ra 3 ) and k (see [19] for how these constants 
can be determined). 

Choosing M — 1, M = 0, g v G {n^+1, — nj, A such that Ajv,jv,i,o = 0, -E = E+C N ^ >lfi , 
we can use the result in Section 14.11 to transform this solution to an eigenfunction of the 
Inozemtsev-type differential operator H^t'~(x,x) with g v — A — g v {v — 0,1,2,3). As 
a suitable integration region we now choose a closed path C in the complex y-plane such 
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that jy l (x, x, y)f(y) is analytic in some neighborhood of this path and such that ( HHj) is 
fulfilled. For example, for p e cjiZ+u^Z, i e {1,2,..., N}, we can choose as C a figure-eight 
contour in the y-pl&ne which encloses y = Xi counterclockwise, y = 2p — Xi clockwise, and 
which does not contain branching points of ^ N fj x (x, x, y) other than y = x,^ 2p—Xi insida^l; 
see Fig. A (the condition fHHj) is satisfied, because the function f(y)-^^I N ^ 10 (x,x,y) — 
^ N 10 (x,x,y)-^f(y) comes back to the original value after analytic continuation along 
CI " 




Fig. A. Figure-eight contour C. 



Obviously there are many more such closed paths C, which we call suitable. We thus 
obtain the following. 

Proposition 4.2. For N , N non-negative integers such that N + N > 0, n u non-negative 
integers (v = 0, 1, 2, 3), let g u e {n u + 1, — n u }, A such that 2XN — 2N = g + g\ + g 2 + g%, 
d v = A/2 — g v , H^'~(x,x) as in [35]) , and ^f^'-(x,x) as in [3d]) . Then, for any constant 
E, there exit constants ti(i — 1, 2, . . . , no + n\ + n-i + n^) and k such that, for any suitable 
closed path C in the complex y-plane (e.g. the figure-eight contour in Fig. A), the function 



$[tl~,(x,x 



nf=i - y)e x {xj+y) 



cU7=i^j-y)^(Pi+y) x m 

(3 \ no+ni+n2+ri3 

\{9 1+v {yf"~ n A J] 9 1 (y + t j )dy, 
u=0 J j=l 

r(+) 



is an eigenfunction of H N ^(x,x) with the eigenvalue E. 



Note that the constants tj and k can be specified by the condition that the function 
( |58|) satisfies ( 1571) with the eigenvalue E = E + C N ^ 10 . Note also that there is another 
expression of solutions of the Heun equation (l57j) by the Hermite-Krichever ansatz [21], and 
we can obtain a similar result to Proposition 14.21 

We now describe a classical example of the Bethe ansatz where the constants tj and k 
can be specified in a simple manner, li n = 1, n x = n 2 = = 0, (|57j) reduces to a Lame 
equation, and its solution in (155]) can be written as 

f(y) = exp (-C(t)y) = exp (r /l t 2 /(2a; 1 )) exp (-&(t)y) (60) 

with t such that p(t) = —E; see [21] and references therein. By specializing to this case and 
<7o = — 1, gi = 92 — 93 — we obtain the following (note that c = if g± = g 2 = g^ = 0). 



2 The second author apologizes for a vague description of the integral routes in [12l Proposition 6]. The 
orientations of them should be specified as the figure-eight contour in this paper. 
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Proposition 4.3. For N a positive integer, N a non-negative integer, let A = (2N + 
1)/(2N), d = A/2 + 1, d x = d 2 = d 3 = A/2, H l +l(x,x) as in $3$), and vfy(x,x) as 

in (36)) . Then, for any closed suitable path C in the complex y-plane as described above 
(e.g. the figure-eight contour in Fig. A ), and for any constant t such that p(t) is finite, the 
function 



*(+) ( ~\ I U.f=M^j - vWijxj + y) My + t ) ; f M \ 

w„f, tea; / — kt exp — d>i (t )y — — — — dy 61 

NJfK } 'cUUhfa-vyoite + v)* d i(y) 2 

is an eigenfunction of H^^(x,x) with the eigenvalue E = —p(t). 
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A Elliptic functions 
A.l Definitions 

The Weierstrass p-function with periods (2cui,2u! 3 ) is defined as follows: 



X' 

(m,n)eZxZ\{(0,0)} 



with VL m<n = 2mwi + 2nu 3 . We also recall the definitions of the corresponding Weierstrass 
zeta- and sigma-f unctions, 

Cte) = -+ ( _o + 77^ + F^M 

X (m,n) G ZxZ\{(0,0)} L % m ' n m ' n m > n > 

and 

a{x) = x II l( 1 "r^) exp (^ + 2^ 

( m ,n)eZxZ\{(0,0)} L V m ' n/ V m ' n m - 

respectively. 

We use the following symbols, 

u = 0, u 2 = -ui - u 3 , (63) 

e u = p(u v ), ■q v = C{u v ) (i/ = l,2,3), (64) 
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g = e , r = — , R = , p = (65) 

LO\ IT 7T1 

where we regard ui as fixed. We also need 6 v (x) = r & u (x/R ) q) [y = 1,2,3,4) with the 
Theta-functions $ u (x,q) as usual [H], i.e., 

oo 

9 1 (x) = 2 ^(-l)« e T " i(n+1/2)2 sin(2n + l)x/R, 

n=0 

oo 

e 2 ( x ) =2J2 e T7ri(n+1/2)2 cos(2n + l)x/R, 

(66) 

9 3 (x) = 1 + 2 eT ™ 2 cos 2nx/R, 

n=l 

oo 

6 A (x) = 1 + 2 ^(-l)™"^™ 2 cos 2nx/R. 

n=l 

We also define 

M*) = jM (* = 1,2,3,4). (67) 
Here and in the following we use the following shorthand notation, 

Q' v {x) = ^-e u ( x ), e v {x) = ^-e v {x) = ?EJLe v ( x ). (68) 

ox Op ojf or 

A. 2 Properties 

We recall some well-known properties of the Weierstrass elliptic functions that we need (see 
e.g. [14J, Chapter XX): 

p(x) = -C(x), C(*) = ^4, (69) 

e 1 + e 2 + e 3 = r] 1 +r]2 + r] 3 = 0, (70) 

p(x + 2u u ) = p(x), ((x + 2u v ) = C(x) + 2 Vu (i/= 1,2,3) (71) 

p(-ar) = p(x), £(-x) = -£(x), a(-x) = -er(a?), (72) 

and 

{((xx) + ((x 2 ) + C(x 3 )) 2 = p(xx) + p(x 2 ) + p(x 3 ) (xi + x 2 + x 3 = 0). (73) 

As will be seen, the last identity plays an important role in the proof of our main result. 
Other identities that we need are the heat equation satisfied by the Theta-functions, i.e., 

6l(x) = 26 u (x) (z/ = 1,2,3,4) (74) 

(obvious from the definitions), and two identities obtained from the well-known duplication 
formula 

0i(2ac) = Ce 1 (x)6 2 (x)6 3 (x)e 4 (x) 



15 



where C is a constant (see e.g. [H], Example 5 at the end of Chapter XXI), i.e., 

01(2*) 1 (e[{x) , 9' 2 {x) + e_^x)_ + 9' 4 (x) 



e 1 (2x) 2 ye^x) e 2 (x) e 3 ( x ) e 4 (x) 

and 

9>{{2x) _ 1 A 6»(x) 1 ^ gjx) 

^(2x) 4^^(x) 2 0„(x) 0„(x) 

Using (1671) and the heat equation we can write these identities as 

4 

2 



I 4 

h(2x) = -J2Mx) (75) 



and 



#i(2a:) 1 Qvi x ) 1 , , ^ , , s , s 

(76) 



We also need the relations between the functions defined in (1671) and the Weierstrass 
zeta function: 

= C(^) - — , 4u+i{x) = t{x + Vu)-Vu-— (^ = 1,2,3) (77) 
(this is a simple consequence of 

B x {x) = C (r(x)exp(--^— ), 0„ + i(a:) = C>(x + £*;„) exp f -^z - J (i/ = 1,2,3) 

for constants C v (y = 0, 1, 2, 3), which can be obtained by comparing the product represen- 
tation of the Weierstrass sigma function in [13], § 20.421, with the product representations 
of the Theta-functions in [H], § 21.3). 



In the following we collect several identities needed in the proof of our main result. 
Proposition A.l. The following holds true, 



and 



for u = 0,1, 2, 3 ; and 



K+i(x) 2 = 2^f( + p(x + u v ) + ^ (79) 



v v +i{x) v^+i{x) U\ 2 

for < n < v < 3, where 

e„,o = e u {y= 1, 2, 3), e 2j i = e 3 , e 3)1 = e 2 , e 3j2 = e 1 . (81) 
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Moreover, 

u=Or=± u=0 U+LK ' 



- y)4>i(x + y) = ~ ^r^24> v+ i(x)4>i(x - ry) 

9\{x — ry) 3r]i 



(82) 



f^ ± O x (x-ry) ux 



and 

/or f = 0, 1, 2, 3 ; and 

^ U>i{x - ry)4>i(x - sz) + ^(y - rx)4>i{y - sz) + 4>i(z - rx)(j>x(z - sy)j 



r,s=± 



2 f Q\{x- ry) 0i(x - rz) 0i{y - rz) 37/i 
^—^\8i(x — ry) Qi(x — rz) 0%(y — rz) 2u x ) 



r=± 



(84) 



Proof. Differentiate ( 1771) in x and use ( 1691) to obtain ( 1781 . 

Differentiate (pTj) in x and use (|67l) again to obtain (f> u (x) 2 = 9"(x)/9 u (x) —<f) r u (x). Insert 
the heat equation in ( 1741) and use ( 1751) to obtain ( 1791) . 

We show (jSPj) . Substitute xi = x + ci>„, x 2 = — x — u;^ and X3 = w M — w„ in (173]) . and 
use C(w M - w„) —rj^-rju and p(w M - u^) = e VjlM to obtain 

(C(x + Lu u ) - ((x + ujfj) + ^ - r] u ) 2 = p(x + uj v ) + p(x + ujfj) + e„ iM . 
Insert into this 

C(x + u v ) - ({x + ujfj) +rj li -rj„ = <p u +i{x) - <f>p+i(x) 
(this follows from (1771) ). expand the square, insert ( T79|) . and obtain an identity equivalent 

to am. 



We show f[8"2"j) . Substitute x\ = x — y, x 2 = x + y, x 3 = — 2x in (173"j) . and use ( 157)) and 
( 177]) to obtain 

(0i (x - j/) + 0i(x + y) - 0i(2x)) 2 = p(x - y) + p(x + y) + p{2x). 
Expanding the square and using ( 1791) this can be written as 



>i x- 



y)(j>x{x + y) = </>i(2x) M x ~ r v) ~ a (0 \ ~ 



9\{2x) \^ 0\{x — ry) 3?7i 



O.i^) ^ (><(■>■ ry) 2u x 
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Insert into this (175]) and 

0i (2a;) _ lJ^6 u+ i(x) 1 ^ (6^ + i{x) 9 v+1 (x) rji e v>l 



e^x) A^e v+l (x) 4 o ^ 3 V0 M+1 (x) e v+1 {x) Ul 2 

E ^+iQe) 3?7i 
n ^ + i(x) 2wi 

i/=0 

(we used (175)1 . (ISO]) and ^ M<i , e„,„ = 2(e x + e 2 + e 3 ) = 0) to obtain fl82|) . 

We show ( 1551) . We substitute xi = x + co u , x% = —y — u u and £3 = — x + y in ( 1751) . Then 

(0„+i(ac) - 0„+i(y) - 0i (a; - y)) 2 = (C(z + w„) - C(l/ + w„) - C(a - 2/)) 2 

= p(x + w„) + p(y + u„) + p(a: - y) 

using (177)) . and it follows from (1791) that 

0„+i( a O0«+i(y) + ^+1 ( x _ y) _ 0v+i(y)0i(a; - y) 

= fl„+i(z) fl„ + i(y) flifo-y) 3?7i_ 

^+1(2;) v+1 (y) Qi{x-y) 2u) X ' 

Similarly we have 

- <p u+1 (x)(j) u+1 {y) + <f> v+1 (x)(j>i(x + y) + (j)u+i(y)(f>i( x + v) 
= 0u+i{x) 9 v+1 (y) 0i(x + y) 3r]i 
9 v+ i{x) u+1 (y) 9 1 {x J r y) 2(V 

by substituting x\ = x + u u , X2 = y + oj u and £3 = — x — y — 2u u in ( 1751) . Sum up the last 
two equalities and obtain ( 1551) . 



We finally show ( 1841) . We substitute x± = x + ry, X2 = —x — sz and X3 = — ry + sz, for 
r, s = ±, in ( 1751) . Then 

(0i (sc + ry) — <pi{x + sz) — 4>i{ry — sz)) 2 = p(x + ry) + p(x + sz) + p(y — rsz) 
and, by using ( J79|) as above, 

4>i{x + ry)4>i(x + sz) + 0i(y + rx)4>±(y — rsz) + <pi(z + sx)(f>i(z — rsy) 

9\(x + ry) Q\ (x + sz) 9\{y — rsz) 3r]i 
9i(x + ry) 6i(x + sz) 9\{y — rsz) 2u\ 

Sum the last equality over r, s — ± to obtain ( 1841) . □ 



B Proof of Source Identity (details) 

We compute the functions defined in ( 1T51) . 

Vj = Y^9u,j(t>v+i{Xj) + ^ ^mjm^A0i(Xj - rX K ), (85) 

1/ K^Jr=± 



and thus 

U 



-J2— d J + ™> W = E — (djVj + vl) (86) 

with dj = d/dXj. Here and in the following, we write Y2 U snor t for X^ =0 , E j short for 
E^=i> etc - 

We compute >V = Wi + W 2 + W 3 with 



J J { v 

+ 2 ^29u,jg^j(pu+i(Xj)(f)^ +1 (Xj) 



the sum of all one-body terms, 



m = E \ E [™>kM'x{Xj - rX K ) + mjm 2 K X 2 MXj - rX K f 

+2 g„ }J m K \(f) v+ i(Xj) 4>x (Xj - rX K )^j 



+2mjm 2 K X 2 MXj - X K )<j> x (Xj + X K ) 



the sum of all two-body terms, and 



w - = E ±{ E (<Mfcj - 1) (»K» + -J + 1) 



2ff 



I/, J 



@v+l(Xj 

8v+i{Xj 

+2 E ^ ^^^) + e^x7) 

Changing summations in the last sum we obtain, after some computations, 

m = T E \—9v,A9v,J - l)p(z + uj v ) + 2g v ,j(\d\ + 2m,A A +1 ™ 

J v \ mj Vv+l{Xj) 

+ yZ\(\d\+ 2mj\)(mj\d\ + 2m 2 \ -l)--mjV eW M e„ iA 



17) 



m = J2Yl E E mjm K m L X 2 MXj - rX K )<^{Xj - sX L ) (89) 

J K^J L^J,K r,s=± 

the sum of all three-body terms. 

To simplify the one-body terms we use f l78l) - fl80l) and obtain 



(90) 
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where we used 

y^ff^j = rnj{\d\ +2mj\), e„ t(t = 0, ^{dp + d v )e Uifl = 

following from (fTBlh (T70|) and (|8T|). 

To compute the two-body terms we use ( 1751) . f[T§j) and (1521) to obtain 

™ = £E £ m KK m J m K^ - 1) - rX K ) + ^ 

+ ^{2g^jm K \ + mjm^A^^+^Xj^pO - rX K ) 

where terms 2^ r =± — tXk) /Oi{Xj — tXk) are cancelled. Symmetrizing the terms in 

the first two lines using [m K + mj)X{mjm K \ — 1) = 7j^, 

2du,fm>K^ + mjm 2 K \ 2 = mjm K \(2d v + (mj + m^)A) 

and 

- rXj) = -r<j) X {Xj - rX K ), 

we can write this as 

m = El ^ k (^ x j - rX ^ + — ) 

J<K r=± [ ^ Wl ' 

+ S ^m J m K X{2d u + (mj + m K )\) (<f) u+1 (Xj) - r<p v+ i(X K ))<j)i(Xj - rX K ) | 
We insert ( 1831) and partly undo the symmetrization, 

J<ifr=± L V 

+ y^ mjm K X(2d u + (mj + m y )A) ^-^T 



.j -rX K 



+ ^2mjm K \(2d u + (rrij + m K )\) — 2ra J m 2 K \ 



J K^J v 



e u+x {Xj) 37?! 
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(a factor 2 in the first term in the last line comes from a r-sum) and, after some computa- 
tions, we obtain 

m = J2J2) -TJKPiXj - rX K ) + mjm K X{2\d\ + A{mj + m^A) ^ 7 '^ 

J<Kr=± I 1\ J - 



K 



Y,Y,^(H-m J )gJ 1 fii^\ (91) 



J v 



+ ( 2 ^jk + 3A (H ~ m j) m j(\ d \ + 2m A) ) 

\J<K J ) 

using the short-hand notation in (fl6l) . 

To compute the three-body terms we symmetrize the summations, 
m= Yl ^iJ^Km L X 2 J2 (MXj - rX K )MXj - sX L ) 

J<K<L r,s=± 

+ MXk - rX L )MX K - sXj) + MXl - rXjMXi - sX K )^j, 
which allows us to use the identity (1511) to obtain 

VI) V* ■> .i^/ te-rA-n) k(Xj-rX L ) h(X K -rX L ) 3r h 

Changing summations again we can write this as 

m = yy y 2 mjmK m L ^y ( d ^- rX ^ + m. 

J K^J L^J,K r=± \ iv J A; \ 

and, by simple computations, we obtain 

^ \Vi{Xj-rX K ) lujrj 

Recalling (EH]) and W = Wi + W 2 + W 3 we add the results in (j90|). (j9l|) and (j92|) and 
obtain 

£ = — \ ~ d " + E - + ^) +EE 7J,^(*J - 

J mj V / J<Kr=± 



with 



S = -<j + 2mj\)(mj\d\ + 2m 2 j\ - 1) + 3A(|m| - mj)mj(\d\ + 2mjX) 

+ y^fajx + Amjm K (\m\ - mj - m x )A 2 ) - + y^mj^ dyd^e^ 
j<k ) Ul J fj,<v 



(93) 
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(a factor 2 in the second term in the second line is from a r-sum). Recalling (jTOj) and (ITT|) 
we can write this as in (TT9l) . Some computations show that Sq in (1931) can be simplified to 
the formula given in (TT5l) . 
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